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m ' 

> ' Introduction 

in 

CN . The objective of this paper is to present a very simple and alternative method 

Tij- ■ to the method of undetermined coefficients (see Birkhoff and Rota [2], Braun 

^ ' [3], Coddington and Levinson [4], Krantz and Simmons [5], Pontryagin [6] and 

Robinson [7]) and also to the annihilator method (see Apostol [1]), both very 
well known, of solving a linear ordinary differential equation with constant real 
^ ! coefficients, P{^)x = /, where / is a real function given by a linear combination of 

■ polynomials, trigonometrical functions and exponential functions products. This 

is accomplished by reducing the given equation to the trivial case in which / is 
a polynomial. Most text books develop the method of undetermined coefficients 
but not the annihilator method, considered to be a bit more sophisticated. In 
some books, it is not proven that the undetermined coefficients strategy is a valid 
one and such strategy is, in some texts, named "the guessing method" or "the 
lucky guess method" . In general, applying either the undetermined coefficients 
or the annihilator method requires a good computations amount. There are three 
advantages in the alternative presented here. First, it reduces the computations 
amount necessary in comparison with either the undetermined coefficients or the 
annihilator method. Second, the proof is quite easy. Third, it exhibits a formula 
to find a particular solution to the given linear ode. 
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Preliminciries 

Given an infinitely differentiable function / : R — ?■ C, then f^^^ indicates the 
A;-th derivative of /. We also write, / = /' = f^^\ f" = f^^\ f" = f^^\ etc. 

Lemma 1. The linear ordinary differential inhomogeneous equation with cons- 
tant real coefficients Oj, < j < n, with at least one Oj, 1 < j < n, not zero, 

a„x(") + a„_ix(""^) + ... + aix' + QqX ^ R , x = x{t) ,teR, 

where the function R = R{t) = hnt"" + ••• + hit + h^, bj G M, < j < n, is a 
non zero real polynomial on the real variable t, has a real polynomial solution 
Q = Q{t) -.R^R. Moreover, 

(a) If oo 7^ 0, then we have degree((5) = degree(i?). 

(b) li k = max{/ '■ aj = , j < I}, then we can admit Q = t'^'^^Qi, where Qi is 
a real polynomial satisfying degree (Qi) = degree (i?). 

Proof. 

(a) Let us solve the pair of equations, 

Q{t) = c„r + ... + cit + co, 
;i) aoQ + aiQ' + ... + a,g(j) + ... + a^g^^^) = R, 

by identifying the coefficient of at the parcels ajQ^^\ j < I, noticing that 
at the others parcels the correspondent coefficient is zero. Fixed the parcel 
ajQ^^\ j < I, a, factor of the coefficient comes up with the computation, 

c„_,+, — = Cr,-i+j{n - l + j){n -l+j- - I + 1)^"', 

and the coefficient looked for is then ajCn-i+j ■ 

Therefore, substituting in (1) the expression Q{t) = Cnt^ + ... + Cit + Cq, and 
its derivatives, we conclude that the coefficient of i""' satisfy the identity, 
where the sum is in decreasing order on j, where j = 1,1 — 1, 0, 

(2) aiCn + ... + ajCn-l+j + ■■■ + «0Cn-i = hn-l , ^ = 0, 1, U 
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By the expressions given by (2) we get the trivial matricial equation, 
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(b) In this case, the equation is a„a;''"'^ + ... + ak+ix''^^^^ = R. By (a), the 
equation anV^'^''^^^^^ + ... + ak+iD = R-, where k + 1 < n, has a solution 
y{t) = Q{t), with degree((5) = degree(i?). Integrating {k + l)-times the 
function y = y(t) and choosing zero for independent term, at each step, we 
reach the desired solution ■ 
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Lemma 2. Let P(^) = a„^ + an-i-§^ + ... + + ao-T, where a^ e R, 
< j < n, and / is the identity operator defined over C°°(M). Let p{t) be the 
characteristic polynomial of Then, if Q = Q{t) G C~(M) and 7 G C, 



/'(|){gwe-} - 

Proof. It is clear that 



^^"^^^^.g(") + ... + + p\^)Q' + p(7)g 



2! 



.7* 



(|-7/){e^*QW} = e^*g'W, 



and iterating we arrive at 



d 



Applying Taylor's formula to pit) about 7 we find. 



fe=0 



A;! 
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Substituting ^ for t and computing the resulting operator at e^Q{t) gives 

fe=0 k=0 

Keeping the above notation and supposing Q and it! real polynomials, and 
7 e C, we have the theorem below. 

Main Result 

Theorem. The equation P(^)a; = R{t)e^* has a particular solution Q{t)e^*, 
with Q a polynomial satisfying the equation, 

(3) ?^!(2)g(») + ... + ?^Mg' + ?^MQ = fl. 

Moreover, depending on 7 we have the following properties, 

(a) If 7 e R, we can suppose Q real and so does the function x{t) — Q{t)e^^. 

(b) If 7 ^ R then Q{t) has complex coefficients and z{t) = Q{t)e^^ is a complex 
solution. If 7 = a + where a and /3 are real numbers, then the functions 
x{t) — Re{^(t)} and y{t) — lm{z{t)} satisfy 

P{j^x = R{t)e'''cos/3t , P{^)y = R{t)e"'sen/3t . 

(c) If ^(7) 7^ then, degree((5) = degree(i?). 

(d) If 7 is a root of multiplicity k then we can choose Q{t) — t^Qi{t), with 
degree((5i) = degree(i?). 

Proof of the Theorem. By Lemma 2, searching for a solution x{t) = Q{t)e'* 
we arrive at 



= R{t)e^ , 



n\ ^ 1! 

and we get equation (3). Hence, applying Lemma 1 and solving an easy diago- 
nalizable linear system we find a polynomial solution, Q — Q{t), of (3). 

Now, we proceed to verify the mentioned properties, 
(a) If 7 is real, it is obvious that we can suppose Q also real. 
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(b) It is also trivial. 

(c) Clearly, if ^(7) 7^ then degree((5) =degree(i?). 

(d) If 7 is a root of multiplicity k then the equation (3) becomes, 

which has a polynomial solution y — Q'^^^ degree((5^'^^) = degree(i?). In- 
tegrating k times the function y{t) but choosing at each one of the k steps 
the independent term as zero, we find a particular solution as desired ■ 

Examples 

(1) Let us solve the ode x'" - bx" + 3a;' + 9x = t^e^* , x = x{t) : R ^ M. 
The characteristic polynomial is 

p(A) = - SA^ + 3A + 9 = (A - 3)2(A + 1) . 
The general solution of the associated homogeneous equation is 

Xh{t) = cie^* + C2te^^ + cse"* , ci, C2, C3 e M . 
By the theorem above, the ode has a particular solution Q{t)e^^ such that, 

S\ ^ ^ 2\ ^ V. 0\ ^ ~ 

Since p' = SX^ - lOA + 3, p" = 6A - 10 and p'" = 6, we easily get 

Q'" + 4Q" = . 

Substituting y — Q" in this last equation we find the ode y' -\-Ay — , which 
has a solution y = ^ + at^ + + ct^ + + e, for certains a, 6, c, d, e e R. 
It is easy to verify that 

_ ht^ _ 15*2 m _ 15 

~ T ~ Te ^ Te ~ ~~u ^ i28 ~ 512 ' 

Hence, by integrating y{t), we can choose 

O- — -— 5t^_15t^ 
^ ~ 168 96 ^ 64 256 ^ 256 1024 ' 
Therefore, the general solution of the given ode is, 

x(t) = cie^*+C2te^*+C3e~*+ f ^ —+- — — i^V^* c's e R. 

W 1^2 ^Vl68 96 64 256 256 1024/ ' ' 
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(2) Let us solve the ode x" + 2x' + 2x = e°^*sin^i , a, ^ e K , x = x{t) : R M . 

The characteristic polynomial is 

p{X) ^ A' + 2A + 2 = (A + 1)2 + 1 - (A + 1 - i)(A + 1 + i) . 

The general solution of the associated homogeneous equation is 

^hif) = Cie~* cost + C2e~* sin t , ci , C2 G M . 

If ^ = then, the given ode is homogeneous. So, the problem is solved. 

If ^ 7^ then, from the identity e"* sin /3t = Im{ 6^"+*'^)* } it follows that the 
imaginary part of a complex solution of the complex ode 

(2.1) x" + 2x' + 2x = where = a + i/3 , 

is a (real) solution of the initial ode. By the theorem proved above, the 
complex ode (2.1) has a solution Q{t)e^* such that Q satisfy the equation, 

(2.2) P-Mq" + pMq' + Mq = Q"+p'{j)Q'+p{j)Q = 1. 
Let us analize three cases. 

(a) Suppose 7 7^ — 1 ± i (hence, 7 is not a characteristic root). Then, the 

constant function Q{t) = l/p{'~f) sstisfy equation (2.2), the complex 

valued function z{t) = p{'y)e"'^ /\p{'y)\'^ solves the complex ode (2.1) and 

the real valued function Im{2;(t)} = Im{p(7)e'*'*}/|p(7)p is a particular 

solution of the given ode. The general solution of the initial ode is , 

/ \ -f -t Im{p(7)e'''* } 
xit) — c\e cost + C2e smt + , , ..^ , ci , C2 G M . 

b(7)r 

(b) Suppose 7 = -1 + i. Wc can write (2.2) as Q" + 2iQ' = 1. Clearly, 
Q' — l/{2i) satisfies such equation. So, we choose Q = —ti/2. Hence, 
the complex valued function z{t) = Q{t)e"'^ = — |e^*ie** solves (2.1) 
and the real valued function Im{2;(i)} = — |e~*cosi is a particular 
solution of the given ode. The general solution of the initial ode is, 

x(t) — cie~^ cost + 626"^ sin t — ^e~*cost, ci , C2 G M . 

(c) Suppose — —1 — i. This case is similar to (b). The solution is 

x{t) — cie~*cosi + C2e~*sini + ^e~*cosi, ci , C2 G IR ■ 
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Introduction 

m 

^ I The objective of this paper is to present a very simple and alternative method 

Cjj . to the method of undetermined coefficients (see Birkhoff and Rota [2], Braun [3], 

■ Coddington and Levinson [4], Krantz and Simmons [5], Pontryagin [7], Robinson 
^ ' [8], Ross [9] and [10], and Zill [11]) and also to the annihilator method (see [10], 

[11] and Apostol [1]), both very well known, of solving a linear ordinary diffe- 
rential equation with constant real coefficients, P(^)x = /, where / is a real 
^ , function given by a linear combination of polynomials, trigonometrical functions 

■ and exponential functions products. This is accomplished by reducing the given 
equation to the trivial case in which / is a polynomial. Most text books develop 
the method of undetermined coefficients but not the annihilator method, consi- 
dered to be a bit more sophisticated. In most books, it is not proven that the 
undetermined coefficients strategy is a valid one and such strategy is, in some 
texts, named "the guessing method" or "the lucky guess method" . It is worth to 
point out that a paper by Ross, see [9], explains very nicely why the undetermined 
coefficients strategy works. In general, applying either the undetermined coeffi- 
cients or the annihilator method requires a good computations amount. There 
are three advantages in the alternative presented here. First, at many times it 
reduces the computations amount necessary in comparison with either the unde- 
termined coefficients or the annihilator method. Second, the proof is quite easy. 
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Third, it shows a very clear algorithm to find a particular solution for the given 
linear differential equation. 

As it is well known, the general solution of the homogeneous equation associ- 
ated to the given ode is a consequence of the Fundamental Theorem of Algebra 
(see Oliveira [6]), and can be found in the textbooks mentioned above. 

Preliminaries 

Given an infinitely differentiable function / : R — > C, then f^'^^ indicates the 
A;-th derivative of /. We also write, / = /(°), /' = f" = f^^\ f" = f^^\ etc. 
Moreover, we indicate C°°(R; C) = {/ : R ^ C : f^^'^ exists ,VA; e N}. 

Lemma 1. Let us consider the linear ordinary differential inhomogeneous equa- 
tion with coefficients Oj € C, < j < n, with at least one aj ^ 0, 1 < j < n, 

QnX^"^ + an-ix^"~^^ + ... + aix' + aox ^ R , x ^ x{t) ,teR, 

where the function R — R{t) — -|- ... -|- bit + bg is a non zero polynomial in 
the real variable t, with coefficients bj G C, < j < n. Then, the equation has a 
polynomial solution Q — Q{t) : R — >■ C. Moreover, 

(a) If gq 7^ 0, then we have degree((5) = degree(i?). 

(b) If A; = max{j : ai — ,1 < j} exists, then we can admit Q — t'^'^^Qi, where 
Qi is a polynomial satisfying degree((5i) = degree(i?). 

(c) If all the coefficients above are real, then Q and Qi are real polynomials. 
Proof. If k exists, then it is clear that k < n. 

(a) Let us solve the pair of equations, 

Q{t) = c„r + ... + cii + co, 
(LI) aoQ + aiQ' + ... + ajQ^^^ + ... + a„Q(") = R, 

by identifying the coefficient of t*^"', / < n, at the parcels ajQ^^\ j < I, 
noticing that at the others parcels the coefficient of t"'~'' is zero. Fixed 
the parcel ajQ^\ j < I, b, factor of the coefficient comes up with the 
computation. 
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(]3 

c„_,+,— = c„_;+,(n - l + j){n -l+j- l)...(n - / + 1)^"', 
and the coefficient looked for is then ajCn-i+j • 

Therefore, substituting in (1.1) the expression Q{t) = c„t" + ... + Cit+Co, and 
its derivatives, we conclude that the coefficient of t"^' satisfy the identity 
(in which, the sum is in decreasing order on j, where j = /, / — 1, 0), 

n! (n — I + j)\ 



[1.2) aiCr, 



+ ... + ttjCn-l+j- 



+ ... + aoCn-i = bn-i , / = 0, 1, n . 
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Prom the expressions given by (1.2) it follows the trivial matricial equation, 
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'■ + Ofe+ii/ = R has a solution y{t) = Q{t), with degree((5) 
= degree(i?). Integrating {k + l)-times the function y = y(t) and choosing, 
at each time, zero for independent term, we reach the desired polynomial. 

(c) It is trivial ■ 

Lemma 2. Let P(^) = On^ + an-i-^Ji + ... + oi^ + oq/, where aj e R, 
< J < n, and / is the identity operator over C°°(M;C). Let p{t) be the 
characteristic polynomial of -P(|). Then, if Q = Q{t) e C°°(R;C) and 7 G C, 
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Proof. It is clear that 

(|-7/){e^*g(t)} = e^*g'(t), 

and iterating we arrive at 

(|-7/)'{e-*g(t)} = e-*g«(t). 

Applying Taylor's formula to p{t) about 7 we find, 

k=0 

Substituting ^ for t and computing the resulting operator at e^^Q{t) gives 

fc=o ■ fc=0 

Let us keep the notation of Lemma 2. Let us consider R, a real non zero 
polynomial, and a pair of numbers (7, 5) such that 

• 7 is complex and 5 is real, with the condition that 5 = if 7 e R . 

Main Result 

Theorem 3. The equation P(^)x = R{t)e^^^^^ has a particular solution Q{t)e"'*^' 
where Q{t) is an arbitrary polynomial satisfying, 

(3.1) l^lt[)Q(») + ... + ?^Mg' + PMg=B. 

Moreover, 

(a) If 7 e M, then we can suppose Q real. So, is real the solution x{t) — Q{t)e^*. 

(b) If 7 ^ M, then z{t) = Q{t)e^^^^^ is a complex solution. If 7 = a + where 
q;,/3 e M, then the functions x{t) = Re[z{t)] and y{t) = Im[z{t)] satisfy 

P{-^)x^ R{t)e'^'cos{/3t + S) , P{^)y ^ R{t)e''' sm{/3t + S) . 

(c) If p(7) 7^ 0, then degree(g) = degree(it!). 

(d) If 7 is a root of multiplicity k, then we can choose Q{t) = t^Qiit), with 
degree (gi) = degree (i?). 
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Proof. Let us apply Lemma 2. Searching for a solution Q{t)e'^'^'^^ we arrive at 



n\ 1! 



and so, wc get equation (3.1). Hence, applying Lemma 1 and solving a trivial 
lower triangular linear system we find a polynomial solution Q{t) of (3.1). 

Now, we proceed to verify the mentioned properties. 

(a) If 7 is real, it is obvious that we can suppose Q also real. 

(b) It is obvious. 

(c) It is obvious. 

(d) If 7 is a root of multiplicity k then the equation (3.1) becomes, 

?q2)Q<~) + ... + 4|2)q(*) = R^t), ^ , 

n! k\ 

which has a polynomial solution y = Q'^^\ <legxee{Q^'^'>) = degree(i?). Inte- 
grating k times the function y{t), but choosing at each time the independent 
term as zero, we find a suitable polynomial solution of (3.1) ■ 

Examples 

(El) Let us solve the ode x" — 2x' -\-2x = t^e* sin(3t + 5), where x = x{t) : M ^ R. 
The characteristic polynomial p(A) satisfies, 

p(A) = - 2A + 2 = (A - 1)^ + 1 , p'{X) = 2(A - 1) and p"(A) = 2 . 

The general solution of the homogeneous associated equation is 

cie* cos t + C2e* sin t , ci , C2 G K. 

By the theorem, the complex equation 

z"-2z' + 2z = 

has a solution z{t) = Q(t)e^^^^'^^*^^\ with Q a polynomial that satisfies 
(El.l) Q"+p'{l + 3i)Q'+p{l + 3i)Q = f . 



Substituting p(l + 3i) — —8 and p'{l + 3i) — 6i in equation (El.l) we get, 

Q" - SQ' + 6iQ = f . 
So, we have that Q{t) = ^ + At + B, with A, B G C It is easy to see that. 



« 2 4. 1 16 



Q(t) = — r --t+ — + —i. 
' 6 9 18 27 



Hence, we find that 

z{t) 



_ if2 _ 4^ , J_ , 16 • 
gfc g*- ^ 18 ^ 27* 



e* [cos(3t + 5) + i sin(3t + 5)] 



lm[z{t)] = e* 



_ cos(3t+5) _ 4t siii(3t+5) , siii(3f+5) , 16 cos(3f +5) 
6 9 18 27 



The general solution of the given ode is, 

x{t) = Cie* cos t + C2e* sint + 



where ci and C2 are arbitrary real constants. 



(E2) Let us solve the ode x'" - 5x" + 3x' + 9x = t^e^\ where x = x{t) : M 
The characteristic polynomial is 

p(A) = - SA^ + 3A + 9 = (A - 3 f{X + 1) . 

The general solution of the associated homogeneous equation is 



cie^* + C2te^* + cae *, ci, 02,036 



By the theorem above, the ode has a particular solution Q{t)e such that, 

^\ ^ ^ 2\ ^ V. 0\ ^ ~ 

Since p' = SA^ - lOA + 3, p" = 6A - 10 and p'" = 6, we easily get 



g'" + 4g" = t^ . 
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Substituting y — Q" in this last equation we find the ode y' -\-Ay — , which 
has a solution y = ^ + at^ + + ct^ + + e, for certains a, 6, c, e e R. 
It is easy to verify that 

5t^ Iht^ Iht 15 

^^^^ ~ T ~ Te ^ Te ~ ~~u ^ 128 ~ 512 ' 

Hence, by integrating y{t), we can choose 

^ 5t=^ 15t^ 



168 96 64 256 256 1024 
Therefore, the general solution of the given ode is, 

x(t) = cie'^'+coie'^ +C3e~*+ \ \ e"** , ds e 

1^2 ^Vl68 96 64 256 256 1024/ ' ' 
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